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Abstract 

An exact inhomogeneous solution of Einstein's field equations is shown 
to be able to inflate in a non-uniform way in the early universe and explain 
anomalies in the WMAP power spectrum data. It is also possible for the 
model to explain the accelerated expansion of the universe by late-time inho- 
mogeneous structure. 
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1 Introduction 

In the following, we shall investigate an exact inhomogeneous cosmological solu- 
tion of Einstein's field equations obtained by Szekeres for an irrotational dust 
dominated universe and subsequently generalized by Szafron [2] and Szafron and 
Wainright [3], H] to the case when the pressure p is non-zero. The cosmological 
model contains the Friedmann-Lemaitre- Robertson- Walker (FLRW) cosmology and 
the spherically symmetric Lemaitre-Tolman-Bondi (LTB) model |5l [6l [7] as special 
solutions. The Szafron model describes a more general inhomogeneous spacetime 
with no prescribed initial symmetry [8]. 

We will show that the Szafron inhomogeneous cosmology in the early universe 
near the Planck time can lead to an inhomogeneous inflationary period, which can 
predict that the primordial power spectrum is not uniform across the sky. Recent 
investigations have revealed that there appears to be a lack of power in the CMB 
power spectrum above ~ 60° and anisotropy in hot and cold spots on the sky P [TUl 
[TTl |T21 [13] • Moreover, there has been the claim that there is a peculiar alignment 
between the quadrupole and octopole moments called the "axis of evil" pll]. These 
anomalies do not agree with the standard homogeneous and isotropic inflationary 
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models [HI [151 IlHl [12] . The anomalies in the WMAP data could be explained by 
the inhomogeneous solution presented in the following. Several possible explanations 
for the anomalies in the WMAP data have been proposed including an exact planar 
solution of Einstein's field equations [18], the violation of Lorentz symmetry and 
rotational invariance [191 (201 [211 [22] and a cutoff of inflation that leads to a loss of 
primordial power spectrum [23] . 

The problem of explaining the acceleration of the universe as determined by 
supernovae data and the cosmic microwave background (CMB) data is one of the 
most significant outstanding problems in modern physics and cosmology [211 [SSI [261 
[27] . The standard explanation is either based on postulating a cosmological constant 
A or assuming that some form of uniform dark energy with negative pressure exists 
in the universe [28]. The explanation as to why the cosmological constant is zero or 
very small has led to a crisis in physics and cosmology. 

The LTB model has been used to explain the non-Gaussian behavior observed 
in the WMAP data [26l [27] and give a possible explanation for the late-time ac- 
celeration of the universe [291 ISOl [SH [32] . A criticism of the LTB model is that it 
assumes a spherically symmetric universe with one spatial degree of inhomogeneity, 
requiring a center of the universe and that observers be located not too far from the 
center to avoid undetected large anisotropy. With this restriction the LTB solution 
has provided a toy-model description of the late-time, non-linear regime with voids 
and collapsing matter. 

2 Inhomogeneous Cosmological Solution 

The metric takes the form0 

ds"^ = dt^ - i?2(x, t){dx^ + dy"^) - ^^(x, t)dz^, (1) 
where -R(x, t) and S'(x, t) are to be determined from Einstein's field equations: 

G^u = Rf,u — -^OfiuR = SnGT^i, + Ag^j^^, (2) 
where T^i, denotes the perfect fluid energy-momentum tensor 

Tf,u = (p + p)Uf,u,y - P9iiu, (3) 

and A is the cosmological constant. Moreover, p denotes the energy- density of 
matter, p the pressure and the velocity field of the fluid which is normalized to 
u^u^ = 1. The coordinates are assumed to be comoving so that = 6q and ii'^ = 
where ii^ = du^/dt. Let us introduce the notation 

i?(x,t)=exp(/3(x,t)), (4) 
^We use units with the speed of light c=l. 
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and 

S{'x,t) = exp(a(x, t)). (5) 

Szafron [2l|3lll] solved the Einstein equations with p and A = 0, generahzing 
the dust solution of Szekeres p]. There are two classes of solution (3' and jS' = 
where (3' = 8(3/ dz. The more general solution /5' 7^ is given for A = by 

i?(x, t) = exp(/3(x, t)) = a{z, t) exp(z/(x)), (6) 

S{yi,t) =exp(a(x, t)) = h{z) exp{—h'{x.))d exp{P{-x,t)/dz, (7) 

where 

exp(-z/(x)) = A{z){x^ + + 2B,{z)x + 2B2{z)y + C{z). (8) 

The A{z)^ Bi{z), B2{z), C{z) and h{z) are arbitrary functions of z. Moreover, a{z, t) 
satisfies the equation 

^ d{z,t) a\z,t) k{z) _ ^ , . , . 

which has the same form as one of the Friedmann equations in FLRW cosmology, 
except that a = a{z,t) and k = k{z). The function k{z) is determined by 

k{z) = Ai^A{z)C{z) - Bl{z) - Bl{z) - ^^). (10) 

The case is singular. Eq.Q can be formally integrated once pit) is specified. 
The density equation is given by 

ci(x, t) + 2/3(x, t) + d2(x, t) + 2/52(x, t) = -47rG(p(x, t) + 3p{t)). (11) 

An algorithm for generating an exact solution is to specify explicitly p = p{t) and 
solve for a{z,t). The metric is now obtained by solving for i?(x, t) and S'(x, t) 
from Eqs.([n]), dZ]) and ([HD- The equation for the density p(x, t) is given by (ITT]) . 

We shall consider, in the following, the simpler Szafron solution with j3' = 
given by 

ii(x.«).exp(/3(x.<)) = ^^^0^. (12) 

5(x, t) = exp(a(x, t)) = X{z, t) + a(t)S(x), (13) 

^ lUiz)ix^ + y^) + Viiz)x + V2iz)y + 2W{z) 
^""^ l + \k{x^ + y^) ■ ^ ^ 

Now /c is a constant, f/(z), Vi(2;), V2(-2) and W{z) are arbitrary functions of z and 
a(t) is determined by the Friedmann equation 

ait) a^it) k _ . s . s 

2-^ + ^ + ^ = -87rGp t . 15 
ait) a'^it) a^it) 
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We can choose without loss of generahty W{z) = and X{z,t) is determined by 

X{z, t)a{t) + \{z, t)a{t) + \{z, t)d{t) = -8nGX{z, t)a{t)p{t) + U{z). (16) 
The matter density equation is given by 



X{z,t) - X{z,t) 



a(t) 



exp(-a(x)) + H\t) + 



k 



SttG 



p(x,t), (17) 



where H{t) = a(t)/a{t). Eqs. ffTSl) and f|T6|) can be solved once the pressure p = p{t) 
is specified, while f lT7|) determines the density p(x, t). The FLRW spacetime is 
obtained when X{z,t) = U{z) = 0. When U{z) = and Vi{z) = V2{z) = 0, the 
model possesses a 3-dimensional symmetry group acting on 2-dimensional orbits. 
The symmetry is spherical, plane or hyperbolic when A; > 0, A; = or A; < 0, 
respectively. 

When k = the line element becomes 

ds"^ = dt^ - a^{t){dx^ + dy'^) - S'^{^,t)dz^. (18) 

A cosmological solution is obtained for a{t) and ^(x, t) when p{t) and p(x, t) are 
specified. 



3 Coarse- Grained Spatial Averaging of Inhomo- 
geneous Model 

For our exact inhomogeneous cosmological model, we are required to carry out a 
volume averaging of physical quantities. We define for a scalar quantity \E' a coarse- 
grained spatial smoothing [331 EH ES] : 

(vi>(x,t))z> = ^ jj^x^-^{^,t), (19) 

where 

Vd = d'x^ (20) 

is the volume of the simply-connected domain, D, in a hypersurface. We can define 
effective scale-factors for our spatially averaged cosmological model: 

Rn{t) = {Ri^,t))n=[^^^Y\ (21) 

Sn{t) = {S{^,t))n=(^^^^y\ (22) 
where Vio is the initial spatial volume. 
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We define the spatially averaged Hubble parameters 



and the effective Hubble expansion parameter 

H,s{t) = ^{2HRn{t) + HsD{t)). (24) 

Consider a congruence of curves with a time-like unit vector with g^'^Vf^Vi, = 1 
and dV^ /ds = 'S/^V^V'^ is the acceleration of the flow lines. Here, and s denote 
the covariant derivative with respect to g^v and the proper time, respectively. The 
metric tensor h^j, is given by 

h^u = S^u + V^V, (25) 

and describes the metric that projects a vector into its components in the subspace 
of the vector tangent space that is orthogonal to V . 
We define the vorticity tensor 

^^lv = y[uV^], (26) 

and the shear tensor ^ 

o'f_iu = Ofj_u — -9h^u, (27) 

where 

V = V{.V^)- (28) 
The covariant derivative of V can be expressed as 

3 ^ ds 
The volume expansion is given by 

e = h^'-e^, = v^v^. (30) 



"^uV^ = 00^, + a^, + - - -j^K- (29) 



The volume averaging of the scalar ^ does not commute with its time evolu- 
tion [331 El [35]: 

(^(x,t))B - 9t(^(x,t))z5 = (^(x,t))^(^(x,t))z5 - (^(x,t)0(x,t))z5. (31) 

For inhomogeneous cosmology, the smoothing due to averaging of the Einstein field 
equations does not commute with the time evolution of the non-linear field equations. 
This leads to extra contributions in the effective, averaged Einstein field equations, 
which do not satisfy the usual energy conditions even though they are satisfied by 
the original energy- momentum tensor. It is the lack of commutativity of the time 
evolution of the expansion of the universe in a local patch inside our Hubble horizon, 
that circumvents the no-go theorem based on the local Raychaudhuri equation 
[37j . namely, that the expansion of the universe cannot accelerate when the weak 
and strong energy conditions: p > and p -|- 3p > are satisfied. 



5 



4 Inhomogeneous Inflationary Model and Microwave 
Background 

Let us now consider the very early universe and for simplicity k = 0. We choose 

p = —Z = const., (32) 
where Z > 0. Then, a de Sitter solution of Eq. ffTSj) is given by 



a{t)=eMHat), Ha = \h^- (33) 



We have from Eqs. fll6l) and ([Tj 

A(z, t) + A(z, t)Ha + A(z, t)Hl = 8nGZ\{z, t) + e^^{-Hat)U{z), (34) 

and 

3 

X{z, t)Hl - X{z, t) + -Hi exp(a(x)) = 47rG'p(x, t) exp(a(x)). (35) 

We now obtain from (fT2|) and ffT3l) the metric 

ds^ = dt^ - exp{2Hat){dx'^ + dy^) - [A(z,t) + exp{Hat)^{x)]'^dz'^ . (36) 

We see that as t ^ oo the universe inflates in the x — y plane but could for S ~ 
be non- inflating along the z direction. For A = and S = 1, we regain the inflating 
de Sitter spacetime metric 

ds"^ = dt^ - exp{2Hat) {dx^ + dy^ + dz^). (37) 

When inflation ends, the metric in the early radiation dominated universe should 
approach a homogeneous and isotropic FLRW metric as t increases with -R(x, t) ~ 
a{t) and the equation of state p{t) = p(t)/3. After decoupling the FLRW spacetime 
is the matter dominated solution with zero pressure, p(t) = 0. 

Let us investigate how the inhomogeneous solution can give rise to a non-uniform 
power spectrum. The power spectrum for homogeneous and isotropic inflation is 
given by 

(5(k),5(q)) = P(A:)53(k-q), (38) 

where 6(k) denotes the primordial density contrast and the translational invariance 
of the inflationary epoch leads to the modes with different wave numbers being 
uncoupled. Let us assume that there exists a vector in the anisotropic direction of 
unit vector n. We assume parity symmetry k — > — k and denote by P{k) the power 
spectrum caused by the anisotropy of primordial spacetime. We have 



P(k) = P(A;) l + /(A;)(k-n)2 , (39) 
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where = |k|, k is the unit vector along the direction of k and we have kept only 
the lowest power of k ■ n. We require that the power spectrum be scale invariant, 
P{k) ~ with f{k) being independent of k, so that f{k) ~ A where ^ is a 

constant. We now have 



P(k) = P(A;) l + ^(k-n 



(40) 



The inhomogeneous spacetime gives rise to correlations between multipole mo- 
ments that normally are zero for homogeneous and isotropic inflation. Deviations 
from homogeneous and isotropic inflation can be parameterized by 



ei 



a{t) - Rpjt) 

m 



(41) 



where a(t) ~ exjp{Hat), Roif) is the smoothed-out scale factor (12T|) and K{t) is the 
averaged cosmic scale: 

(42) 



K{t) = -{2Rn{t) + a{t)). 



5 Acceleration of the Late-Time Matter Domi- 
nated Universe 

Let us now consider the matter dominated solution with zero pressure, p = 0. We 
have 

a{t) a^{t) k 



+ 



0. 



a(t) a?{t) a^{t) 
This equation has the standard Friedmann solution for k = 0: 

The matter equation is given by flTTI) where 

exp(— «(x)) 

Let us define 



(43) 



(44) 



A(^,t) + (-j S(x) . 



?>H\t) ' 



(45) 



(46) 



and for a spatially fiat universe with = 0: 



^^x(x,t) 



it) 



exp(— a(x)) 



(47) 
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We have 

n^{x,t) + nxix,t) = i. (48) 

We define spatially averaged flm and fix- 

VL„,D{t) = (fi„(x,t))B = / d^Xy/^Q^{^,t), (49) 

nxoit) = {nxi^,t))n = ^ f d'xy/^nxi^,t), (50) 
It follows from (jlHl) that 

nrnDit)+nxDit) = l. (51) 

We obtain when A = 0, the standard FLRW density relation for an Einstein-de 
Sitter universe: 

nn,{t) = 1, (52) 

where 

^m{t) = ^P{t). (53) 

By choosing the parameterization: 

n% = 0.28, n% = 0.72, (54) 

where fl^D ^^"^ ^xd denote the present matter and inhomogeneity densities, we 
can fit the WMAP data l27j and the supernovae data 

6 Conclusions 

We have derived a primordial inhomogeneous inflationary model from an exact inho- 
mogeneous solution of Einstein's field equations. Because the model has a preferred 
direction of anisotropic inflation, it is possible to explain detected anomalies in the 
WMAP power spectrum such as the apparent alignment of the CMB multipoles 
on very large scales and the loss of power for 6 > 60°. The observed CMB tem- 
perature anisotropies can give a window on the primordial inflationary era. We 
have given generic predictions expected from the existence of a preferred inhomo- 
geneous direction during inflation. The inhomogeneous inflation should reveal itself 
in a scale-invariant manner, and lead to predictions for the primordial fluctuation 
correlations determined by a single amplitude A and a unit vector n signifying a 
preferred anisotropic direction on the sky. 

A possible explanation of the accelerated expansion of the universe at late times 
and the WMAP data is attributed to the observed late-time non-linear, inhomo- 
geneous galaxy and void structure. Our exact inhomogeneous solution avoids the 
anti-Copernican assumption of a center to the universe, and it also avoids the pos- 
tulates of an ill-understood cosmological constant and "dark energy" . 
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